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1 .  Introduction 

Since  Landau's  phenomenological  theory  of  the  superfluids,  it  is  well 

known  that  liquid  helium  at  low  frequencies,  long  wavelengths  and  low 

2 

temperatures  can  be  described  by  Landau's  two- fluid  equations.  These 

equations  involve  not  only  the  normal  and  superfluid  density  and  thermodynamic 

functions,  but  also  the  various  sounds  and  transport  kinetic  coefficients  in 

liquid  ^He  and  ^He .  Concerning  the  kinetic  coefficients,  Landau  and 

Khalatnikov^  first  investigated  the  thermal  conductivity  and  viscosity,  which 

were  also  studied  later  by  others.^  Diffusive  processes^  in  liquid  ^He  and 

He-  He  mixtures  have  been  extensively  studied;  Reuppeiner  et  al  measured 

the  thermal  conductivity  (K)  ,  thermal  diffusion  ratio  (K,^)  and  mass 

diffusivity  (D)  as  a  function  of  He  number  concentration,  and  Behringer  and 

Meyer^  investigated  diffusive  relaxation  processes  in  the  normal  and 

superfluid  phase  in  liquid  ^He-^He  mixtures.  All  these  quantities  mentioned 

above  are  closely  dependent  on  the  elementary  excitation,  which  is  appropriate 

3  4 

to  describe  the  given  Bose  system.  Motivated  by  these  studies  of  He-  He 
mixtures ,  in  the  present  paper  we  evaluate  the  temperature  dependence  of  the 
diffusion  constant  (D)  and  thermal  diffusive  ration  between  the  phonon 

and  roton  gases  in  consideration  of  the  scattering  of  phonons  by  rotons  in 
two-  and  three-dimensional  liquid  He. 

Recently  we  have  obtained  microscopically  the  Landau- type  excitation 
spectrum,  which  is  (anomalous)  phonon-like  at  low  momentum  and  roton- like  for 

g 

large  momenta  in  the  ring  diagram  approximation.  Using  this  excitation 

9  10 

spectrum  we  have  successfully  analyzed  the  sounds,  sound  attenuations, 

11  4 

kinetic  coefficients,  etc.,  in  two-  and  three-dimensional  liquid  He. 

Through  the  calculations  we  can  consider  the  excitation  spectra  in  both 


dimensions : 
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£(p)  -  CqP[1  +  7?^  -  +  ...1  .  (2D)  (1.1) 

£(p)  -  CqP[1  +  -  52?^  +  ...1  .  (3D)  (1.2) 


and 


(P-Pq)^ 

£(p)  -  d  +  - ,  (2D  and  3D) 


(1.3) 


where  P  is  the  roton  momentum,  is  the  sound  velocity,  and  A,  P^  and  are 
the  roton  parameters  (energy  gap,  momentum  and  effective  mass,  respectively). 
The  subscript  zero  refers  to  T  -  0  K.  All  coefficients  in  Eqs .  (1.1)  and 
(1.2)  are  positive,  which  are  determined  by  the  potential  parameters.  Here, 
we  have  adopted  a  soft  potential  with  a  Lennard-Jones  type  tail,  which  helps 
to  make  a  smooth  connection  of  the  attractive  part  with  the  soft  repulsive 
core.  In  the  next  section  we  evaluate  the  differential  phonon-roton 
scattering  cross  section,  and  calculations  of  the  diffusion  constant  and 
thermal  diffusion  ratio  are  described  in  Sec.  3.  Finally,  the  results  and 
discussion  are  presented  in  Sec.  4. 


2.  Phonon-roton  differential  scattering  cross  section 

To  obtain  the  two-dimensional  phonon-roton  scattering  cross  section,  in 

this  section  we  consider  a  roton  in  the  presence  of  the  phonon  field.  We  can 

A 

treat  this  roton  as  a  particle  in  a  moving  liquid  He.  Then  there  appears  an 


<  - 

O  . 


additional  term  -P»v,  which  can  be  written  in  symmetric  form 


1 

•“(P*v  +  vP) 
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□ 

□ 


(2.1) 


Codes 


iOist 


- - .....,/or 

Speolal 
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where  P  and  v  are  the  momentum  and  velocity  operators.  The  phonon  field 
changes  the  density  of  the  medium,  and  thus  we  may  expand  the  roton  energy 
[Eq.  (1.3)]  in  terms  of  the  density  p'  to  second  order  (p'  -  P-Pq)  as 


H 

r 


“ro  f; 


21  ,  2  “  * 
dp 


(2.2) 


where  is  given  by  Eq.  (1.3).  Since  the  value  of  most  rotons  is  close  to 

Pq,  we  may  neglect  (P-P^)  and  replace  P  by  P^.  We  may  also  drop  the  term 

(^  p'),  which  is  much  smaller  than  (2.1).  Then  the  interaction  energy 
dp 

between  phonons  and  rotons  can  be  written  as 


V 


i 


(P*v  +  v*P)  + 


1  ^^0  2 


p' 


(2.3) 


where  the  terms  in  the  second  bracket  of  Eq.  (2.3)  have  magnitudes  on  the 

-1  13 

order  of  10  to  1  in  three-dimensional  liquid  helium. 

When  the  roton  changes  momentum  P  to  P'  ,  it  absorbs  a  phonon  with 

^  * 

momentum  p  and  emits  a  phonon  with  momentum  p'  .  Here  we  may  consider  two 
intermediate  processes;  (I)  P  +  p-+P'  -  P  +  p  -  p'  and 

(II)  P  -  p'  -►?'  -  P  -  p'  +  p.  Since  the  roton  momentum  is  much  larger  than 
that  of  the  phonon,  we  may  view  this  interaction  as  similar  to  that  of  heavy 
and  light  particles.  The  momentum  and  energy  conservation  law  in  collision 
processes  can  be  written  as 

'P  -  ‘P'  (|p  +  p  -  p'l  ■  ■ 


(2.4) 
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2 

Under  the  conditions  p  .P'  «  Pq  and  E  —  cp  «  3  /ic  ,  Eq.  (2.4)  becomes 

p  -  P'  -  — ~ — r  {P„m*(ph  -  p'fi'))^ 

2mcPo' 

p2  2 

-  Im-(h-n'))^  .  (2.5) 

where  m,  n  and  n'  are  the  unit  vectors  directed  along  Pq  ,  p  and  p', 
respectively.  Therefore  energy  conservation  implies  p  =  p' .  This  means  that 
particles  do  not  change  the  magni'ude  of  momentum,  but  change  its  direction. 
Taking  account  of  p,p'  «  Pq  and  P  -  Pq  ,  the  matrix  element  given  in  the 
differential  cross  section 


da 


(^)l<F|H3ll>|^  S(Ep-Ei) 


dP 

(2»rH)^ 


can  be  written  in  second- order  perturbation  as 


(2.6) 


<F|H3|I> 


{  [  (m*h)  +  (m*h')  ] 


(n»h')  +  ~(m»h)^(m»n' )^ 

HC 


+  A) 


(2.7) 


A 


2  .2.  ,  ap- 


(2.8) 


Substituting  Eqs.  (2-7)-(2.8)  into  Eq.  (2.6)  and  performing  the  integration 
over  p' ,  we  obtain 
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Aa  - 


i>2  3 

^oP 


l[(m.n)  +  (m»n'))(m«n')  +  — (m*n)^(m»n' +  A^}d5  .  (2.9) 

j  z  z  uc 


SwM  PqC 


Averaging  Eq.  (2.9)  over  all  directions  of  roton  momenCum,  we  finally  obtain 


Aa  - 


^qP  2  1  Pf)  2  4  4 

,  -  ((1  +  cosV>)cos  ^  +  7Tr( — )  (35cos  +  3sin  V 

°  ^  Pq^ 


P  A 

2  2  1  0  2  2  2 
+  30cos  Vi  sin  ^)  +  ~  ^J^(3cos  0  +  sin  V*)  A  )  d^fr  , 


where  is  the  angle  between  the  incident  and  scattered  phonons . 

Through  the  similar  calculations  we  can  obtain  the  differential 

4  14 

scattering  cross  section  of  phonons  by  rotons  in  bulk  liquid  He: 


Aa 


PqP  22  21  ^0  2  23 

(  ;  )  {T(1  +  CQsi))cos^\p  +  r“( — )  (1  +  8cos  ^6  +  -cosV>) 


+  t4(— )(1  +  2cos^^i)  +  A^)  dn 

iZ  /iC 

where  dfl  is  the  solid  angle. 


(2.11) 


3.  Diffusion  constant  and  thermal  diffusion  ratio 

4 

Since  the  excitation  spectrum  in  liquid  He  consists  of  phonons  and 
rotons,  we  may  treat  liquid  He  as  a  mixed  phonon-roton  gas.  In  the 
temperature  range  from  -0.7  to  -1.2  K,  where  the  scattering  of  phonons  by 
rotons  is  dominant,  we  consider  the  diffusion  of  phonons  and  rotons  and  also 
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the  thermal  diffusion  ratio  under  the  existence  of  a  temperature  gradient. 
There  thermal  processes  in  the  mixture  of  a  phonon-roton  gas  are  very  similar 
to  those  of  light  particles  in  heavy  particles.  Since  the  roton  momentum  is 
much  larger  than  that  of  the  phonon,  we  may  treat  the  interactions  between 
phonons  and  rotons  as  similar  to  those  between  light  and  heavy  particles. 
Therefore,  we  may  assume  that  when  a  collision  occurs  only  between  a  phonon 
and  roton,  the  roton  is  at  rest  and  the  phonon  momentum  changes  only  its 
direction  with  the  same  magnitude. 

Let  n  and  n  be  the  number  densities  of  the  phonon  and  roton, 
p  r 

respectively,  and  d  the  angle  between  the  direction  of  a  phonon  with  momentum 

p  and  the  x-axis.  Equation  (2.11)  represent  the  differential  scattering  cross 

section  that  the  momentxjun  p  of  a  phonon  will  be  changed  to  the  momentum  p' 

directed  into  the  solid  angle  element  dO,  The  collision  probability  of  a 

phonon  per  unit  time  becomes  n^cda,  where  c  is  the  phonon  velocity.  Let  the 

number  of  phonons  in  a  given  unit  voliane  with  momentum  between  p  and  p  +  dp 

3 

and  directed  into  the  solid  angle  dfl  be  n(p,5,x)d  p.  Then  the  total  number  of 
phonons  changing  p  to  p'  in  dO  and  conversely  p'  to  p  in  dfl'  are, 
respectively,  given  by 

d^p  J  n^cn(p,«,x)  da  , 

d^p  I  n^cn(p' ' ,x)  da 

The  difference  between  the  above  two  processes , 

d^p  n^c  J[n(p',5',x)  -  n(p,«,x)]  da  ,  (3.3) 


(3.1) 

(3.2) 
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represents  the  change  in  the  number  of  phonons  in  the  volume  element  d  p, 
which  must  be  equal  to  the  total  time  derivative  of  n: 


'^phr^">  - 


c  cos^ 


3x 


-  n  c 
r 


|[n 


(c' ,9' ,x)  -  n(p,e,x)]  do' 


(3.4) 


If  the  temperature  and  concentration  gradients  vary  very  slowly,  we  can 
express  the  phonon  distribution  function  as 


n  -  n^Cp.x)  +  5n(p,«,x) 


(3.5) 


where  Sn  is  a  small  correction  term  to  the  equilibrium  distribution  function 
n^.  Now  we  assume  that  5n  is  linear  in  the  gradient  of  the  concentration  and 
temperature,  and  has  the  form 


6n  “  costff(p,x) 


(3.6) 


Substitution  of  Eq .  (3.6)  and  cosd'  -  costfcos^  +  s  in^sin^cos  (0-0' )  in  Eq. 

(3.4)  and  then  integration  over  the  solid  angle  yields 


4  2  4  4 

a  PokgT^ 

■'phr^">  ■  2  6  P 

4irl<  PqC 


(3.7) 


P  PA 

r  -  ^  ^(^)2  ^  ^(^) 

45  25>c''  9>c  ^ 


+  a: 
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Here  we  have  made  use  of  the  relation  p  -  ak  T/c,  where  a  is  a  constant^^  (see 

B 

the  Appendix).  From  Eqs .  (3.4)  and  (3.7),  we  obtain  the  following  expression 
for  f (p,x) : 


f(p,x)  -  - 


/  2  6 
1.  >’0'=  io 

"t 


(3.8) 


The  diffusion  flux  i  along  the  x-axis  is 


-1 


3  r  2  3 

i  -  I  nc  COS&  d  p  -  c  cos  0  f(p,x)d  p 


(3.9) 


Substituting  Eq.  (3.8)  into  Eq.  (3.9)  and  integrating  over  momentum  p  gives 


i  - 


/  2  7  . 

Un\i  PgC  ^ 

uV.rn  v.y 

0  r  B 


(3.10) 


Let  the  phonon  concentration  of  the  mixture  be  rj  -  np/(nQ+n^)  -  n^/n.  Then 
the  diffusion  flux  i  becomes 


i  - 


/  2  7 

^0  fill  ^  d_,  .  .31, 


(3.11) 


Comparing  Eq.  (3.11)  with  the  expression^^  for  the  diffusion  flux 


i  -  -nD(7f,  +  ^  VT)  . 


(3.12) 
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we  obtain 


/  u  2  7 

4jf>(pQC 

3a*p2rn_.k‘;‘ 


-  rjT  f^ln(n) 


(3.13) 


(3.14) 


Since  the  number  densities  of  phonons  and  rotons  per  unit  volume  can  be  given 
by 


n 

P 


2!r(3)  4,r( 


_v 

^■kY'C 


(3D) 


(3.15) 


n 


r 


2P^(ak 

(2jrK^)^/^ 


-VkgT 

e 


(3D) 


(3.16) 


we  obtain  tne  diffusion  coefficient  and  thermal  diffusion  ratio  of  phonons  in 
bulk  liquid  helium: 


(27r)^/%^(Xc)^  q.,  W 

-  4  4  1/r  9/2  ^  « 


(3.17) 


•5  ^  (3).  -1,  ,1/2  „2  3,1  ^  T^-5/2 

3  +  (  yf^)  (’T#*)  PqC  (2  +  ® 

B 

[1  . 


(3.18) 


Through  similar  calculations  as  done  above,  we  obtain  the  diffusion 
coefficient  and  thermal  diffusion  ratio  in  thin  helium  films  as 


11 


D(T)  - 


^-7/2  //V 


(3.19) 


K(T) 


2  +  {r(2)l‘^(27r^)^/^PQC^(^  +  ^)(kgT) 

B 

[1  +  {f(2))'^  (2jrM)^/^PQC^(kgT)‘^/^e 


-3/2 

e 


-vv^2 


(3.20) 


where  we  used  the  number  densities  of  phonons  and  rotons  per  unit  surface 
given  by 


n 

P 


k  T 

LL21  ,  V  2 

2n 


P«  n  iiy  -A/k_T 

n  -  -2t  (2^mV)^/2  ^  B 


(3.21) 


(3.22) 


and  r'  is  given  as 


r'  -  J  +1^  (-7)^ 

4  32  ^^c, 


MC 


+  2A 


4.  Results  and  discussion 

Since  the  elementary  excitation  spectrum  in  liquid  helium  consists  of 
the  combined  phonon- roton  gas,  we  may  consider  the  diffusion  of  phonons  into 
the  roton  gas,  and  also  may  include  the  thermal  diffusion  under  the  existence 
of  the  temperature  gradient.  We  may  treat  the  diffusion  and  thermal  diffusion 
of  the  phonons  in  the  roton  gas  as  similar  to  those  of  light  particles  in 
heavy  particle...  For  convenience  we  neglected  phonon-phonon  or  roton-roton 
scatterings,  permitting  only  phonon-roton  scatterings. 
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The  thermal  diffusion  ratio  represents  the  ratio  of  the  thermal 

diffusion  coefficient  due  to  the  temperature  gradient  to  the  diffusion 

coefficient  given  by  the  concentration  gradient.  The  diffusion  coefficient 

must  be  positive,  but  the  diffusion  ratio  may  be  positive  or  negative.  The 

negative  value  of  the  coefficient  means  that  the  heavy  particles  will 

appear  in  the  low- temperature  region  and  the  light  particles  in  the  high- 

3  4  17 

temperature  region.  In  He-  He  mixtures  the  ratio  K,^  becomes  negative  near 
the  plait  point  and  positive  near  transition  temperature.  Since  thermal 
diffusion  ratios  in  our  calculations  are  positive,  phonons  and  rotons  are 
located  in  the  low-  and  high-temperture  regions,  respectively. 

The  diffusion  constant  and  thermal  diffusion  ratios  of  the  phonon  gas  in 
bulk  liquid  helium  and  in  thin  helium  films  are  given  by  Eqs.  (  3  .  17 )  -  <  3 . 20)  . 

The  potential  and  roton  parameters  are  chosen  by  the  analysis  of  the 

18  8 
excitation  spectnam  for  the  bulk  case  and  the  specific  data  for  the  thin 

helium  films  (Table  I).  For  the  density  variations  of  the  roton  parameters, 

14 

we  have  adopted  the  results  obtained  by  Khalatnikov,  who  analyzed  the  data 
19 

of  Dietrich  et  al  given  by 


-1.142, 


2  =2. 
dp 


-2.06, 


(f~)  (■ 
0 


dp' 


1 

3 


The  above  parameters  in  thin  helium  films  are  not  known  and  thus  are  assumed 
to  have  the  same  values  in  the  bulk  case.  The  constants  a  and  a'  in  Eqs. 
(3.7)  and  (3.19)  can  be  determined  by  the  calculation  of  the  collision  time 
for  the  phonon-roton  scatterings,  which  characterize  the  thermal  conductivity 
and  viscosity.  The  numerical  values  of  these  a  and  a'  are  6.28  and  4.67, 
respectively. 
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Figure  1  illustrates  the  diffusion  constant  D  as  a  function  of 

temperature  in  the  bulk  and  for  films.  As  temperature  increases,  the  roton 

density  increases  more  rapidly  than  the  phonon  density.  Therefore,  the  mean 

free  path  of  phonons  becomes  shorter  and  the  coefficient  decreases  rapidly. 

The  temperature  variations  of  the  coefficients  K.^,  are  shown  in  Fig.  2. 

Below  -  0.6  K  the  coefficient  is  almost  constant.  In  the  temperature  range 

from  -  0.6  to  -  1.1  K  the  coefficient  ratio  increases,  which  means  that 

phonon-roton  scatterings  play  a  main  role  in  the  transport  phenomena.  This 

temperature  range  agrees  well  with  the  results  obtained  by  Landau  and 
12 

Khalatnikov.  For  temperature  far  below  -  0.6  K,  the  phonon  population  is 
dominant,  and  thus  phonon-phonon  interactions  becomes  important,  while  roton- 
roton  interactions  and  five-phonon  processes  will  be  dominant  above  -  1.1  K. 

In  thin  helium  films  the  temperature  range  in  which  phonon-roton 
scatterings  are  effective  and  dominant  becomes  -  0.3K<T<0.8K.  This 
range  is  in  good  agreement  with  the  results  obtained  from  the  evaluation  of 

9 

kinetic  coefficients  in  thin  helium  films.  All  arguments  raetioned  with 
respect  to  the  bulk  helium  will  be  valid. 
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Making  use  of  Eqs.  (3 . 13) - (3 . 15) ,  (3.21)  and  (All)  from  Ref.  15,  we  can 
easily  obtain  the  characteristic  time  ^  expressed  as  a  function  of  the 
phonon  momentum  and  temperature,  respectively,  as 


-1 

’^ph-r 


^0^  1  9  ^0  2  ^0  2 

N  — [/  +  T7  (~)  +  —  A  +  2A^] 


(Al) 


and 


-1 

^ph-r 


61061 
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2  3  3 


9.  IQ  2 
32 


flC 


A  + 


2 
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(A2) 


where  is  the  number  of  rotons  per  unit  area.  Setting  (Al)  equal  to  (A2)  , 
we  get 


pc  -  [ 


61061,1/3 

3!r(3)J 


V 


(A3) 


Therefore,  we  obtain  a'  -  4.67  in  two-dimensional  liquid  helium. 

Through  a  similar  calculation  we  can  evaluate  the  characteristic  time 
expressed  by  temperature  and  phonon  momentum: 


-1 

^ph-r 


-  4wc  N 
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4x)f(  pc 
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tic 


+  A^] 
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and 
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-1 

’^ph-r 
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P  PA 

A  1  0  2  2  0  2 

"  f  (;:r)  ^ 


(A5) 


From  Eqs.  (A4)-(A5)  we  obtain  a  -  2jr  in  the  bulk  case. 
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Table  I .  Theoretical  Parameters 

T,  A/kg  (K)  ^ 


3D 

2.18  X  10'^  a'^ 

8.616 

1.930 

0.153mHe 

238.21 

2D 

2.79  X  10'^  A'^ 

4.12 

1.02 

0.77mHe 

164 . 4 
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Figure  Captions 

1.  Diffusion  constants  of  phonons  in  the  roton  gas  as  a  function  of 

temperature.  The  solid  and  dotted  lines  represent  the  theoretical  curves 

4 

in  two-  and  three  dimensional  liquid  He. 

2.  Thermal  diffusion  ratio  of  phonons  in  the  roton  gas  versus  temperature. 

The  solid  and  dotted  lines  represent  the  theoretical  curves  in  two-  and 

4 

three-dimensional  liquid  He. 
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